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1.  Introduction 


Cross-linked  polymer  networks  swollen  with  polymeric  solvent  have  shown  adaptive  mechanical 
response  (/).  This  frequency-dependent  response  makes  these  gels  desirable  for  applications 
such  as  tissue  simulants  for  ballistic  testing. 

This  experimental  observation  is  qualitatively  modeled  with  the  discrete  slip-link  model  (DSM) 
introduced  by  Schieber  and  Schieber  et  al.  (2,  3)  and  further  developed  by  Khaliullin  and 
Schieber  (4-6).  Two  applications  of  the  DSM  are  used.  First,  an  ideal  entangled  network  (IEN), 
described  by  Jensen  et  al.  (7),  is  modeled.  This  network  is  stoichiometrically  cross-linked  and  no 
dangling  ends,  soluble  sub-structures,  or  solvent  chains  are  present.  The  entanglements  are 
permanent  (or  trapped),  since  all  chains  are  attached  to  the  permanent  network  at  the  ends. 

Despite  being  trapped,  relaxation  through  monomer  shifting  can  occur  between  entanglements. 

Secondly,  the  DSM  is  applied  to  these  ideal  entangled  networks  swollen  with  solvent  chains. 

This  application  is  a  coupling  of  two  architectures:  (1)  entangled  network  chains  (ENCs)  and  (2) 
entangled  solvent  chains  (ESCs).  The  experimental  system  modeled  is  a  stoichiometrically 
cross-linked  polydimethylsiloxane  (PDMS)  network  in  the  presence  of  a  high-molecular  weight 
PDMS  linear  solvent.  The  network  relaxation  and  flow  behavior  are  predicted. 


2.  The  Discrete  Slip-Link  Model 


2.1  Ideal  Entangled  Network 

As  illustrated  by  references  7  and  8,  the  network  chain  can  be  described  and  coarse-grained  by 
random  walk  statistics  and  modeled  as  a  single-chain  mean-field  model.  Entanglements  define 
the  primitive  path  of  the  chain.  The  parameters  and  variables  of  the  coarse-grained  chain  are 
listed  in  table  1.  In  this  application,  Z,  chosen  from  a  distribution  for  a  given  chain,  is  constant 
since  entanglements  are  trapped  in  the  IEN  model.  Monomer  density,  TV*,  fluctuates  due  to  Kuhn 
step  shuffling  through  the  slip-links  in  the  sliding  dynamics  (SD)  process. 
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Table  1.  Parameters  and  variables  of  the  coarse-grained  chain. 


Description 

Parameters 

K 

Total  number  of  Kuhn  steps* 

(3 

Entanglement  spacing  parameter 

tk 

Average  time  of  transferring  one  Kuhn  step  through  an  entanglement 

Variables 

Z 

Number  of  entangled  strands 

Nt 

Number  of  Kuhn  steps  in  strand  i 

Qi 

Orientation  vector  between  entanglements  i  and  i  +  1 

rCD 

i 

Characteristic  lifetime  of  i  entanglement  determined  by  constraint  dynamics  (CD) 

*7  =  N,  S  indicates  that  we  can  work  with  two  different  molecular  weights  in  the  second  model  application. 


The  random  walk  of  Nt  Kuhn  steps  in  a  strand  between  two  entanglements  with  a  connector 
vector,  Qi ,  has  a  Gaussian  free  energy  given  by 

FsiQuNj)  _  3 Q\  3 

kBT  2Nia%r  2 

where  kB  is  the  Boltzmann  constant  and  T  is  the  temperature.  The  free  energy  of  the  network 
chain  (NC)  composed  of  Z  strands  is  a  sum  of  the  free  energies  associated  with  the  entangled 
strands: 

z 

FKC{tt)  =  Y,F8{Qi,Ni),  (2) 

i= 1 

where  Q  =  (Z,  { ArJ ,  {Qi},  { TpD }  j  is  the  cross-linked  chain  conformation.  The  equilibrium 
distribution  of  such  a  chain  is  given  by  the  modified  Maxwell-Boltzmann  relation  (5) 


2n  NiCi^ 


(1) 


piqm  = 


^Zk  Ef.i  ffi. 

j 


exp 


r  F(Q)  i 

rv(z-i)i 

kBT  \ 

exp 

kBT  \ 

z-\ 


IIpcd(Td: 


(3) 


i= 1 


where  j/  r)  (  r^'r)  j  is  the  probability  density  for  the  ilh  entanglement  with  a  characteristic  CD 
lifetime,  TpD,  determined  by  CD.  It  is  determined  self-consistently  from  the  destruction  of 
entanglements  in  a  system  by  SD.  Entanglements  between  two  ENCs  will  never  be  destroyed  by 
CD;  hence,  t{'V)  is  infinitely  large  and  for  those  pCD(rCD)  =  2d(l/rCD)/(rCD)2.  The  Kronecker 
delta  function  5(i,j )  :=  t)t.3  is  responsible  for  the  conservation  of  Kuhn  steps  in  a  chain  and 
J  =  (1  +  l//?)^-1  is  the  normalization  constant.  The  entanglement  chemical  potential  of  the 
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surrounding  chains,  he,  is  responsible  for  fluctuations  in  the  number  of  entanglements  on  the 
chain. 

2.2  Entangled  Network  Chains  and  Entangled  Solvent  Chains 

We  combine  two  architectures  of  chains:  one  for  ENCs  and  one  for  ESCs.  The  solvent  chains 
will  have  two  free-ends  where  entanglements  can  be  created  and  destroyed  by  SD,  while  some 
fraction  of  entanglements  can  be  created  and  destroyed  anywhere  on  the  strand  by  CD.  This  latter 
process  is  caused  by  the  SD  of  surrounding  solvent  chains.  While  still  another  fraction  of 
entanglements  cannot  be  created  or  destroyed  by  CD,  as  these  are  entangled  with  the  network 
chains  and  they  are  trapped.  The  ENC  will  have  a  fraction  of  entanglements  created  and 
destroyed  by  CD  due  to  entanglements  with  ESC. 

2.3  Stress  Relaxation  of  Network  Chains  and  Solvent  Chains 


The  stress  tensor  is  obtained  from  the  chain  free  energy 


r7(t) 


fdF{  ft) 

\^Q~ 


T,{Ni},{QjfH} 


7 


(4) 


where  (...)  is  an  ensemble  average.  The  relaxation  modulus,  Gy(t),  is  obtained  using  the 
Green-Kubo  expression 

=  n  {rxy{0)TXy(t))eclr/  ,  (5) 

The  subscript  7  is  either  "NC"  or  "SC"  for  solvent  chain.  The  relaxation  modulus  for  each 
structure,  Gnc  and  G sc  is  obtained  from  Green-Kubo  simulations  of  the  two  ensembles,  while 
G(t)  for  the  combined  system  is  a  linear  superposition  of  the  relaxation  moduli  of  the  two 
components 

G(t)  =  w  ncG  nc  (t)  +  wscGsc(t),  (6) 

where  wnc  and  wsc  is  the  volume  fraction  of  network  chains  and  solvent  chains,  respectively. 
The  network  chains  and  solvent  chains  influence  each  other  through  CD. 


2.4  Equilibrium  Dynamics  for  ENC  and  ESC 


At  equilibrium,  Kuhn  step  shuffling  between  entangled  strands  is  due  to  Brownian  forces  and  free 
energy  differences.  The  evolution  equation  for  the  probability  density  describing  the 
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conformation  of  a  chain  is  given  by  the  Chapman-Kolmogorov  equation  (2) 


dp(Q,t  |  n0,f0) 
dt 


-£ 
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k  •  |  n0,fo) 


jk(q  |  n')P(n',t  \  n0,t0)dn' 


w(tf  |  n)P(n,t  \  no,tQ)drt,  a) 


where  k,  (Vv)^  is  the  transpose  of  the  velocity  gradient.  We  assume  affine  entanglement  and 
cross-link  motion,  dQjdt  =  n  •  Qr  The  probability  density  of  having  conformation  Q  at  time  t 
given  that  at  t0  the  conformation  was  Qq  is  p(D,  t  \  The  function  W (O  |  fT)  is  the 

transition  probability  from  conformation  to  conformation  Q'  per  unit  time  and  has  five 
contributions 


W 


W$\  +  W|D  +  llfD  +  WdCD  +  W?D, 

1=1  J 


(8) 


where  VFs*h  is  the  transition  rate  probability  for  Kuhn- step  shuffling  through  entanglement  i.  kkf D 
and  ITp30  represent  destruction  and  creation  of  the  entanglements  on  the  end  of  the  solvent  chain 
due  to  SD  (these  terms  are  zero  in  the  IEN  and  ENC  models),  and  W^D  and  W^D  represent 
destruction  and  creation  of  the  entanglements  anywhere  on  the  strand  due  to  CD  (these  last  two 
contributions  are  zero  in  the  IEN  model).  When  deriving  expressions  for  the  transition 
probability,  we  satisfy  detailed  balance 


w{vt'  |  n) 


w(n  |  n') 


Peq(tt) 

Peq(^) 


(9) 


Details  of  these  transition  probabilities  can  be  found  in  Jensen  et  al.  (7).  It  is  assumed  that 
cross-links  are  fixed  in  space. 


2.5  Constraint  Dynamics  Spectrum 

It  is  necessary  to  obtain  a  CD  spectrum  pCD  (r)  such  that,  by  assuming  binary  entanglements,  the 
destruction  of  entanglements  by  SD  and  CD  will  be  equal  in  time.  This  can  be  done  by  keeping 
track  of  the  entanglements  from  their  individual  time  of  creation  to  their  individual  time  of 
destruction.  By  doing  this,  we  weight  short-lived  entanglements  more  than  long-lived  ones. 
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This  gives  us  the  cumulative  distribution  of  lifetimes,  P(f),  related  to  p(  r> 


through 


(T 


CD 


) 


fd(t) 


r 

Jo 


m 

pCD(rCD )  ( 

tcd  exp  (- 


r)  drCB 


r 

Jo 


r 

Jo 


Jo 

P§??(r0D) 

TCD 


)drCD 

exp  (-^n)  drCD 


r 

Jo 


pCD(TCD) 

— run — 


drCD 


(10) 


because  /0°°  6(l/r)/r3dT  =  0. 


Since  we  later  compare  the  theory  with  experimental  PDMS  data,  we  perform  simulations  with 
Nk  =  307  and  (3  =  25.3.  The  estimate  of  these  parameters  is  described  in  section  3.  The 
resulting  fd(t)  from  a  solvent  chain  simulation  without  CD  is  shown  in  figure  1. 


t/r*- 


Figure  1.  fd(t)  obtained  from  a  simulation  of  solvent  chains  with 
A^k  =  307  and  j3  =  25.3  and  an  ensemble  of  100  chains. 
The  crosses  are  simulated  data,  while  the  full  line  is  a 
fitted  spectrum:  to  =  1,t\  =  5000,  T2  =  1000000, 
oi  =  0.5,  and  a 2  =  0.3. 
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The  following  spectrum  is  used  to  fit  well  the  CD  normalized  lifetime  distribution  probability  of 
an  ideal  entangled  network  with  solvent  chains: 


PCD(t)  =  wsc 


Tai  1H{t  —  Tq)H(ti  —  T )  r"2  1H(t  —  T\)H (t~2  ~  t)t{ 


CX.\  — OL2 


+  WNC 


ot-1 

m/r) 


+ 


a  2  <*2 

r2  —tx  ^ a\—a2 
/  1 

Ot-2  1 


+ 


a  i 


cz  2  oc  2 

r2  —t1  “  ^ai—a2 
/  1 

Ot-2  1 


T 


(ii) 


The  delta  function  assigns  an  infinite  lifetime  to  a  fraction  jonc  of  entanglements. 


3.  Results 


3.1  BSW  Spectrum 

To  estimate  the  Fourier  transform  of  G(t),  we  fit  the  simulated  data  with  a  mathematical 
expression  given  by 


G(t)  =  (G%-G0)  jf  ^exp^  )dr  +  G0 


(12) 


where  G%  is  the  entanglement  plateau  and  G0  is  the  equilibrium  plateau  and  h{r)  a  modified 
Baumgaertel,  Schausberger,  and  Winter  (BSW)  spectrum  (9) 


h{r)  = 


-  t)H{t  -  Ti)  nj"i  T) 


m  Tf.  Z—Ti- 1  l  ii  1  Qj— aj+i 


Era 

i= 1 


nr-  v 


(13) 


The  fit  is  performed  by  minimizing  the  y2  function  assuming  constant  relative  error 


^2 _ ^  —  g(Xi) 

i= 0 


(14) 


where  /*  is  the  numerical  value  obtained  at  x%  point  and  g(xi)  is  the  analytical  function  that  is 
used  to  fit  the  data.  For  the  IEN  results,  we  fit  only  one  mode,  and  equation  5  reduces  to 

KT)  =  T_  aH(T 1  “  T)H(r  ~  To)  (I5) 

T1  ro 

where  T\  and  r0  are  the  longest  and  shortest  relaxation  times,  respectively. 
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To  transform  between  the  experimental  time  domain  to  the  frequency  regime,  a  one-sided  Fourier 
transform  is  used,  G*(u>)  :  =  iuT  {(?(£)} 

3.2  DSM  Parameters 

To  calculate  the  relaxation  modulus,  experimental  parameters  are  necessary  including  the 
molecular  weight  of  the  polymer  and  Kuhn  length,  Mw  and  MK,  respectively;  polymer  density,  p; 
plateau  modulus,  G°N;  and  temperature,  T.  The  data  for  the  system  are  presented  in  table  2.  The 
polydispersity  index,  MW/MN,  is  also  given  as  a  measure  of  the  quality  of  the  melt. 


Table  2.  Experimental  data  for  PDMS. 


Name 

Mw  (kg/mol) 

Mk  (kg/mol) 

Mw  /Mn 

p  (g/cm3) 

T  (K) 

G%  (kPa) 

T314 

314.3 

0.381 

1.15 

0.965 

293 

227 

Reference 

(?) 

(?) 

(?) 

(?) 

(?) 

(?) 

The  two  input  parameters  for  the  DSM  model  are  the  number  of  Kuhn  segments,  TVk,  and  the 
model  parameter,  /3.  NK  can  be  calculated  using  ArK  =  MW/MK.  The  data  in  Ressia  et  al.  (10) 
are  used  to  find  j3  by  fitting  the  Fourier  transform  of  the  following  expression  to  it: 


G(t)  —  G°n  ^  exp  dr 

where  G°N  is  taken  from  Khaliullin  and  Schieber  (5) 

0  pRT  ( (  (3  Y*-1  nk-P-2\ 

N  mw  ^U  +  iy  /5  +  i  J 


(16) 


(17) 


The  results  are  shown  in  figure  2. 
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Figure  2.  T314  blue  symbols  are  experimental  data  from  Ressia  et 
al.  {10)  fitted  with  a  two  mode  BSW  spectrum. 


With  these  parameters  and  following  the  procedure  outlined  in  Chantawansri  et  al.  (S),  the 
relaxation  modulus  is  transformed  to  the  frequency  domain  and  compared  with  data  as  seen  in 
figure  3.  tk  was  also  retrieved  when  this  fit  to  the  experimental  data  was  done. 

The  values  for  the  two  parameters  are  shown  in  table  3 


Figure  3.  T314  blue  symbols  are  experimental  data  and  the  line  is 
the  DSM  prediction  with  tk  =  0.002/rs  . 
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Table  3.  Parameters  needed  for 
DSM  program. 


Name 

Nk 

P 

tk 

T314 

825 

25.3 

0.002 ps 

3.3  Ideal  Entangled  Network 

In  the  ideal  network,  we  have  only  trapped  entanglements. 

In  figure  4,  we  show  the  Green-Kubo  results  of  simulating  100  IEN  chains  for  Nk  =  307  and 
(3  =  25.3.  The  results  are  normalized  with  the  entanglement  plateau  modulus,  which  is 
determined  analytically  as  (2,  3,  5) 


G°N  _  /y\  _  +  (3 

pRT/M  ~  '  ;eq  “  [3+1 


(18) 


Figure  4.  Relaxation  modulus  produce  by  the  IEN  model  with 

A^k  =  307  and  (3  =  25.3.  The  red  line  is  the  BSW  fit  of 
the  DSM  prediction. 


The  resulting  G*  spectrum,  for  Nk  =  307  and  6  =  25.3,  is  shown  in  figure  5.  The  difference  in 
magnitude  of  G'  and  G"  makes  it  difficult  to  see  the  shape  of  G[  which  is  rendered  in  the  blow  up 
in  the  figure  inset,  revealing  the  three  frequency  regimes.  The  BSW  parameters  are  listed  in  table 
4. 
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Figure  5.  Dynamic  relaxation  spectrum  for  Arj<  =  307  and 

fj  =  25.3  obtained  with  the  IEN  model.  G'  is  several 
orders  of  magnitude  larger  than  G"\  the  shape  of  G'  can 
be  seen  in  the  inset. 


As  discussed  by  Jensen  et  al.  (7),  experimental  observations  of  stoichiometrically  imbalanced 
networks  show  that  G"  and  G1  are  of  the  same  order  of  magnitude  at  intermediated  frequencies; 
the  DSM  suggests  that  energy  dissipation  is  largely  a  result  of  the  dangling  ends  and  soluble 
structures. 


Table  4.  BSW  parameters  obtained  by  fitting 
equation  6  combined  with  equation  4 
to  the  DSM  prediction  presented  in 
figure  2. 


Nk 

P 

Co 

a 

TO 

n 

307 

25.3 

7.5426 

0.0141 

0.5 

2706 

3.4  Entangled  Network  with  Solvent  Chains 

The  DSM  is  used  to  qualitatively  capture  the  trends  of  experimental  data  at  times  scales  slower 
than  glassy  modes.  For  a  quantitative  description,  better  characterization  of  network  structure 
and  solvent  polydispersity  is  required.  An  example  of  experimental  data  is  shown  in  figure  6  (7). 
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Figure  6.  G'  for  a  PDMS  gel  containing  50%  solvent  of 
molecular  weight  139,000  g/mol. 


We  simulate  network  chains  and  solvent  chains  with  the  same  parameters,  i.e., 

TVk  =  307,  /3  =  25.3,  and  iusc  =  0.1.  The  Green-Kubo  results  from  these  simulations  are  shown 
in  figure  7  together  with  the  corresponding  Gien  obtained  in  the  previous  section.  The  length  of 
time  for  relaxation  increases  in  the  presence  of  solvent  chains,  as  expected.  When  compared, 
Gnc  clearly  shows  the  CD  effect  in  comparison  to  Gien  where  G0  is  reached  much  quicker  and  is 
increased.  The  presence  of  solvent  chains  increases  G°N.  The  analytic  expression  for  G°N  is 


pRT/Mk 


-JjW  (ZENc)eq  +  -jys  ((^ESc)eq  “  1 
WNC  +  3  Wsc  f  +  3  \ 

N£  3+  1  +  N*  V  P  +  1  J 


(19) 


Mk  is  the  molecular  weight  of  a  Kuhn  step  while  (ZENc)eq  and  (ZEsc)eq  are  the  average  number 
of  entangled  chains  for  the  network  and  solvent  chains,  respectively,  while  and  are  the 
total  number  of  Kuhn  steps  in  the  network  chains  and  solvent  chains,  respectively.  The  top 
horizontal  line  in  figure  7  is  an  analytic  estimate  of  the  plateau  modulus,  Gq,nc  given  by 


Go,nc 

pRT/M 


0.74Zt  +  1 


(20) 


Z'i  is  the  average  number  of  trapped  entanglements: 


Zt  —  iunc  ((-^enc)  eq  —  1) . 


(21) 
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The  bottom  horizontal  line  in  figure  7  represents  the  estimated  analytic  plateau  modulus,  GoAna, 
of  the  network  with  solvent  chains,  given  by 


Go, ana  _  __  Go,NC 

pRT/M  ~  WNC  pRT/M 


wnc(0-74Zt  +  1). 


These  later  two  equations  suggest  the  G0  is  a  function  of  trapped  entanglements  only. 


(22) 
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Figure  7.  Green-Kubo  simulation  results. 

iVN  =  jyS  =  307,  p  =  25.3,  and  wsc  =  0.1. 


G(t)  for  the  blend  is  fit  again  as  described  previously  and  converted  into  the  frequency  domain  as 
shown  in  figure  8. 

Figure  8  shows  the  G*  result,  and  it  is  seen  that  G'  and  G"  are  becoming  closer  in  order  of 
magnitude  compared  to  the  ideal  entangled  network  with  no  entangled  solvent  present.  A 
qualitative  comparison  of  this  result  with  experimental  data  from  Mrozek  et  al.  (7)  is  shown  in 
figure  9. 

A  number  of  assumptions  listed  in  Jensen  et  al.  (7)  related  to  requiring  better  characterization  of 
the  experimental  network  and  solvent  chains,  suggest  an  explanation  regarding  the  discrepancy 
between  the  DSM  prediction  and  the  data.  Additionally,  a  comparison  with  experimental 
polydisperse  solvent,  which  is  used  for  results  in  the  next  section,  and  the  DSM  prediction  of  the 
same  monodisperse  solvent  is  shown  in  figure  10,  making  apparent  an  effect  of  polydispersity  on 
mechanical  properties. 
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Figure  8.  G*  result  for  N ^  =  N ^  =  307,  /3  =  25.3,  and 
u>sc  =  o.l. 


Figure  9.  DSM  prediction  for  PDMS,  N ^  =  N^_  =  307,  (3  =  25.3 
and  ruse  =  0.1  and  experimental  data  for  a  network  (1 17 
kDa,  PDI=1.55)  with  10%  polydisperse  139  kDa  solvent. 


3.5  Simple  Elongation  of  Ideal  Entangled  Network  with  Solvent  Chains 

The  DSM  is  used  to  predict  the  mechanical  response  of  an  ideal  entangled  network  swollen  with 
monodisperse  entangled  solvent.  For  a  simple  pull  test  for  constant  elongation  rate  of  5  in/sec 
with  an  initial  sample  length  of  1  in,  the  elongation  ratio,  A  (A  =  l/L),  is  given  by 

A  =  At  +  1  (23) 
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Figure  10.  DSM  prediction  for  monodisperse  solvent, 

N =  492,  j3  =  25.3,  and  PDMS  experimental  data, 
187.5  kDa,  PDI=1.6. 


where  A  =  5  1/s.  For  the  DSM,  the  elongation  ratio  is  related  to  the  logarithmic  (Hencky)  strain, 
which  is  used  by  the  algorithm,  by 

e  =  log  (A)  (24) 

and  consequently  the  true  strain  rate  can  be  expressed  as 


.  _  A  1 

Xt  +  1  t  +  l/X 


(25) 


For  the  melt  component  of  the  system,  the  longest  relaxation  time  is,  Td  =  0.03  s,  obtained  from 
Sliozberg  et  al.  (12).  The  resulting  Weissenberg  number,  Wi  =  erd,  is  pictured  in  figure  11. 
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Figure  1 1 .  Weissenberg  number  for  uniaxial  extension  for  PDMS 
melt,  187.5  kDa,  A  =  5  1/s,  Td  =  0.03  s. 


This  result  is  consistent  with  the  DSM  prediction,  shown  in  figure  12,  which  shows  that  the 
solvent  chains  are  not  experiencing  nonlinear  deformation.  For  this  system,  the  main 
contribution  to  the  stress  is  from  the  network  component.  The  system  response  is  a  weighted 
average  of  the  network  and  solvent  stress,  depicted  by  the  blend  prediction.  To  further  illustrate 
this  minimal  contribution  to  stress  from  the  melt,  an  additional  result  is  plotted  in  figure  12  for  a 
prediction  of  just  50%  network  solution  diluted  with  non-entangled  ideal  solvent,  for  which 
j3 solution  _  5o.6  was  calculated  using  a  scaling  argument  regarding  the  entanglement  spacing  for  a 
polymer  chain  in  solution  according  to 

omelt 

^ solution  i_ _  (26) 

0 

where  0  is  the  volume  fraction  of  the  network.  It  can  be  seen  that  the  two  results  are  very  nearly 
the  same. 
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Figure  12.  True  stress,  a,  vs.  stretch  ratio,  A  for  uniaxial  tension 
(A  =  5  1/s)  using  the  DSM  for 

=  307,  iV|  =  492,  (3  =  25.3,  and  wsc  =  0.5.  For 
the  diluted  network,  fjsolvlum  =  50.6,  =  307. 


A  qualitative  comparison  with  experimental  data  from  Sliozberg  et  al.  (72)  is  shown  in  figure  13. 
Again,  characterization  of  network  architecture  and  solvent  polydispersity  is  necessary  for  a 
quantitative  comparison.  Both  of  these  latter  features  could  work  to  bring  down  the  prediction  of 
stress  approximately  0.57  times  making  the  agreement  better  as  can  be  seen  in  the  figure. 

In  a  Mooney  plot,  shown  in  figure  14,  for  the  DSM  predicted  ideal  system  and  the  experimental 
system,  a  linear  mechanical  system  response  can  be  seen. 


4.  Conclusions 


Two  applications  of  the  discrete  slip-link  model  have  been  presented.  The  first  is  used  to  model 
ideal  entangled  networks.  The  other  application  of  the  model  is  a  combination  of  ENCs  and 
ESCs.  As  a  starting  point,  equilibrium  and  flow  experimental  data  were  described  qualitatively, 
and  for  simplicity,  network  chains  and  solvent  chains  were  assumed  monodisperse,  linear,  and  of 
the  same  molecular  weight.  Additional  detailed  characterization  of  the  network  and  melt  is 
necessary  to  make  quantitative  predictions. 
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Figure  13.  True  stress,  a,  vs.  stretch  ratio,  A  for  uniaxial  tension 
(A  =  5  1/s)  using  the  DSM  for 
N £  =  307,  JVf  =  492,  (3  =  25.3,  and  wSc  =  0.5. 
Experimental  data  from  Sliozberg  et  al.  (12)  for  A  =  5 
1/s,  containing  50%  T204. 


A 

Figure  14.  Mooney  plot  for  experimental  data  from  Sliozberg  et 
al.  (12)  for  A  =  5  1/s,  containing  50%  T204  and  DSM 
prediction  for  =  307,  iVj|  =  492,  (3  =  25.3,  and 
ruse  =  0.5  for  uniaxial  tension  (A  =  5  1/s). 
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List  of  Symbols,  Abbreviations,  and  Acroynms 


BSW 

Baumgaertel,  Schausberger,  and  Winter 

CD 

constraint  dynamics 

DSM 

discrete  slip-link  model 

ENC 

entangled  network  chain 

ESC 

entangled  solvent  chain 

IEN 

ideal  entangled  network 

NC 

network  chain 

PDMS 

polydimethylsiloxane 

SC 

solvent  chain 

SD 

sliding  dynamics 
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